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Abstract 

The three integrable two-dimensional Henon-Heiles systems and their integrable per- 
turbations are revisited. A family of new integrable perturbations is found, and TV- 
dimensional completely integrable generalizations of all these systems are constructed by 
making use of si (2, K) © /13 as their underlying Poisson symmetry algebra. In general, the 
procedure here introduced can be applied in order to obtain ./V-dimensional integrable 
generalizations of any 2D integrable potential of the form V(q\,q2), and the formalism 
gives the explicit form of all the integrals of the motion. Further applications of this 
algebraic approach in different contexts are suggested. 
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1 Introduction 



The Henon-Heiles Hamiltonian 



H = 




) 



(1.1) 



was introduced in [1] in order to model a Newtonian axially-symmetric galactic system, and 
it was soon considered as the paradigm of a system that exhibits chaotic behaviour (see for 
instance [HEld]). Later, when its following generalization containing adjustable parameters 
was considered, 



it was proven that the only Liouville-integrable members of this family of generalized 2D 
Henon-Heiles Hamiltonians were given by the three following choices of the parameters 
(see [5] _ [in] and references therein): 

• (i) The Sawada-Kotera case: /3 = 1/3 and $7 = 0. 

• (ii) The KdV case: (3 = 2 and Vt arbitrary. 

• (iii) The Kaup-Kupershmidt case: /3 = 16/3 and $1 = 15 5. 



Here we have used the terminology of [12], since these three integrable cases correspond 

precisely to the stationary flows of three integrable fifth-order polynomial nonlinear evolution 
equations. We also recall that in the KdV case the integral of the motion is quadratic in 
the momenta, while in the Sawada-Kotera and Kaup-Kupershmidt cases is a quartic one. 
Nevertheless, the separability of the latter Hamiltonian has been also been demonstrated 
through an involved canonical transformation in [17\ [T8] , and the explicit integration of the 
three systems has been also studied (see the recent paper [19] and references therein) . There- 
fore, Henon-Heiles (hereafter HH) systems are privileged benchmarks for the study of the 
transition between integrable and non-integrable regimes from different viewpoints, includ- 
ing also their quantizations [2j[3lll]. As a consequence, higher dimensional generalizations 
of the HH systems have been considered [15} \19\ l20"j [2~T] as well as several integrable pertur- 
bations of these Hamiltonians have also been found (see, for instance [22, 23J as two recent 
representative works for the KdV case). 

The aim of this paper is to present a novel Poisson-algebraic approach to the integrability 
properties of the HH systems and their perturbations, that can be straightforwardly applied 
to any 2D integrable potential of the form V{q\,q2)- The cornerstone of this construction 
is obtained by realizing that such 2D Hamiltonians can always be defined as a specific N = 
2 symplectic realization of an 'abstract' Hamitonian % defined on the direct-sum Poisson 
algebra s/(2,R) h%. Afterwards, if we consider a suitable TV-dimensional (herafter ND) 
symplectic realization of 7~L, we obtain an ND completely integrable system whose integrals 
of the motion will be explicitly known. This construction can be interpreted as a Poisson- 
algebraic generalization for the procedure given in [201 [2^] and is based on the more generic 
framework of the so-called coalgebra symmetry approach to integrable Hamiltonian systems 
(see |25]-|29|). In particular, the sZ(2,R) coalgebra symmetry has been recently exploited 
in [30] in order to introduce (super)integrable Kepler-Coulomb and oscillator potentials on 
several ND curved spaces. 

The structure of the paper is as follows. In the next Section the three cases of integrable 
HH systems and all their known integrable generalizations are reviewed. Note that in the 



^ (2) = g (Pi + VD + 6 li + + q 2 2 + a {qjq 2 + (3 q 3 2 ) 



(1.2) 




2 



multiparameter framework given by (|1.2p the oscillator term can also be considered as a 
perturbation, since the oscillator constants 8 and f2 are explicitly taken into account as 
independent parameters. We also present a new perturbation (to the best of our knowledge) 
for the KdV case, that can be obtained as a suitable superposition of previously known 
results. In Section 3 the sZ(2,R) © /13 Poisson algebra symmetry of all the 2D HH potentials 
is shown. Based on this result, in Section 4 we describe the algebraic procedure in order to 
obtain ND integrable generalizations of the HH systems, and we show how the sZ(2,R) © h% 
symmetry will give us automatically the corresponding integrals of the motion. In Section 
5 the explicit expressions for the generalized ND HH systems are given, showing that new 
additional centrifugal/monopole terms coming from the most generic symplectic realization 
of sZ(2, R) © /13 can be included without breaking the complete integrability of the system. 
The final Section is devoted to comment on several possible applications and generalizations 
of the algebraic approach here presented, including the construction of quantum and curved 
analogues of the HH systems. We also stress that the method here presented can be easily 
used to get integrable ND generalizations of any 2D integrable potential of the form V(qf, (72), 
and the explicit example of the Holt potential is sketched. 



2 Integrable perturbations of 2D Henon-Heiles systems 

In this Section we review all the known integrable perturbations of the 2D integrable HH 
systems. We will write all the cases under the same type of parameter conventions in order 
to provide a unified approach to the known results, which are scattered in the literature under 
different notations. 



2.1 The Sawada-Kotera case 



The only known perturbation of the HH-Sawada-Kotera Hamiltonian is obtained by adding 
a centrifugal term to the original unperturbed system, namely |31j 



-H^ = \ij>l+pl) + 5{ q l + q l) 



(2.1) 



and the integral of the motion (which is quartic in the momenta) reads 

Z (2) 



1 2 2 , e/r 2 2, \ 1 2a5q^ , 2 „ 3\ 



+a I aq\q\ 



92 Ji 
2 3 



a 
18 



q\ +piqiq2 



P2<?2 



2pi<?i 



q 2 

+ y [ 2 Pi<72 + PiqiP2 j 



.pi 4a 



(2.2) 



Note that we have ordered the terms in the integral with respect to the perturbation pa- 
rameters. This convention will be preserved from now on, and will facilitate the immediate 
obtention of the corresponding integrals when some parameter(s) vanish. In particular, the 
integral of the motion for the unperturbed system is obtained by taking the A — > limit. 
Note also that the isotropic oscillator term can be easily removed by taking 5 — > 0. 
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2.2 The Kaup-Kupershmidt case 

In this case the only integrable perturbation is given by two rational terms, namely 

^ (2) = \ (Pi + Pi) + %? + 16d) + a f g 2 , 2 + ^ +4 + 4 ( 2 - 3 ) 
i \ 6 j q 1 q 1 



whose constant of the motion is 



j(2) 



l pi + 6(ql[35q 2 1 +p 2 1 ]+2 P 2 1 qj) 

+a (q 2 (q 2 p 2 ~P2Pm) ~ aqj 
3 



6 + <?2 



+ 2q 2 (p\ql-5qi 



+A 



9? 



+ 2ag 2 + — 



2ag 2 + 25 + 4 
9i 



P? + 



2A v 
1 + qJ 



(2.4) 



which is also quartic in the momenta. We recall that the q~ 2 perturbation was given in [20] 
while the term q~ 6 was considered in [32l [331 [M] . Again, the four parameters can be supressed 
independently and the corresponding integral of the motion would be readily obtained. 



2.3 The KdV case 



This is the integrable HH system that admits the largest family of integrable perturbations, 
probably due to its connections with stationary and travelling wave reductions of KdV equa- 
tions [T21 [22l [23l EU [35] . Moreover, under all the perturbationes here considered the integral 
of the motion remains to be quadratic in the momenta. 



2.3.1 The q 2 perturbation 

As usual, a centrifugal term can be added also in this case: 

■H {2) =\{p\+ pi) + Sq 2 + (5 + n) q 2 + a [q 2 q 2 + 2 ql) + A ( 2 .5) 
Z q 1 

and the integral is 

? (2) = S (^p 2 + (35 - n)q 2 ^) - ^p 2 + a {-q 2 p 2 + aq 2 (| + gg) + p 2Pl q^ 

+2a5q 2 qj + A (35 _ ft _ 2aq 2 ) . (2.6) 
Qi 

Note that when all the perturbation parameters vanish, it suffices to take one of the momenta 
as the integral of the motion. This system was considered in |31j and is also a particular case 
of the system II in [20] by taking g$ = g± = = = and g 2 = 5, 4g 2 = 5 + £1, gs = a/4. 
Note that the anisotropy in the oscillator term can be independently removed by taking 

n -> o. 

2.3.2 The Ramani series of polynomial deformations 

The Ramani potentials Vm{qi,Q2) are homogeneous polynomial potentials with degree M 
and given by [331 EE] 

V * = E 2i " 2fc (* fc )rt k l2- 2k i = l,2,... (2.7) 
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namely, the first potentials of this infinite family are given by 



V = l 


(2.8) 


Vi = 2q 2 


(2.9) 


V 2 = Aq 2 2 + q\ 


(2.10) 


V 3 = 8q% + \q\q 2 


(2.11) 


V 4 = IGqi + \2q\q\ + qj 


(2.12) 


V 5 = 32g| + 32qfql + 6qfq 2 . 


(2.13) 


= 2 (Pi + ^i) + "2V 2 + a 3 V 3 


(2.14) 



By considering 



with a 2 = 6, a 3 = a/4 we get the KdV-HH Hamiltonian with Q = 35 y A = 0. But it is well 
known that the full Ramani series can be superposed by preserving the complete integrability 
of the system (see \33\ [36] ) and therefore we obtain the integrable Hamiltonian given by 

M 



u { $ = ^(pi+p!)+X> v * 

1=1 

1 M ^ /• JA 

1 / 2 i 2\ i „ ni-2k 1 K \ 2k i-2k , n i r\ 

= 2 (?i +P2 ) + l^l^ ai2 ( k J 91 9a ( 2 - 15 ) 



i=l fc=0 

whose integral of the motion is: 

M 

-<■-- - _„,,,- .__ „. ,,. .... ,,- » V/ _j 



z m = -12P1 + qwm + Qi^2 aiVi 

i=i 

= -^ + ^ + ^£5^ (2.16) 
2.3.3 Rational perturbations 

In the particular case 5 = Q = 0, a = 1/2 (therefore, with no oscillator potential) the following 
integrable perturbation of rational type is known |23j 



2i+2 

t=0 g i 



\(p\+vD + +^+Et 23i ~ 2k M\ k )J^ (2 - 17) 

i=0 k=0 v 7 9i 



where Vi are the potentials in the Ramani series. The generic integral of the motion for this 
system is given by 

4 2) = QmP2- q 2 P i+ 1 - q hl+ 1 - q t + cql-ql^22^^ 

1 1 R 2 /,' i 1 _A „i+l-2fc 



2i+4-2/c ' 

i=0 fc=0 x ' q i 



(2.18) 
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In particular, when R = 2 we get the integrable system 

^2 = «(Pl +P2) + ^QiQ2 + q 2 +2cq 2 + — + 8fi— j- +32f 2 g (2.19) 

z yz y #1 ?i 9i 

and for R = 3 we recover the system (II) in [20] provided that #i = c, 52 = 0, 53 = 1/8, 2£o = 
\,h 3 = 8£i,/»4 = 326- 



2.3.4 New integrable perturbations 

But it turns out that all the two latter perturbations of the KdV case can be superposed. 
This novel result can be immediately checked by direct computation, and stresses again the 
ductility of the KdV-HH system, that can be perturbed in different directions by preserving 
its integrability. In particular, the rational perturbations can also be included in the full 
Ramani series provided that M > R. This encompasses all the previous KdV perturbations 
with Q = 35 and leads to the new system: 



M R 
& 1/1 . . . . £7i 



M [|] 



1 q i i=ife=o V fc / 

+ EE^" 2fc r , V&32T (M>H). (2.20) 

The corresponding integral of the motion can also be obtained by direct computation and 
reads 



t(2) 2 , 2A 1 2 I -v> Vj+i 

^M.ii = + 91P1P2 J ^2 + Ql 2^ ~ 2^ 7i '~ 



9i s «i 



2 2A 
+ giPiP2 — 2 ^2 

M t^ 1 



\ i=l fc=0 ^ 7 

^ Iff A" 4- 1 fr\ „ l + 1 -2fc \ 

„-_i 1,-n \ / ?i / 



i=l fc=0 



As a particular example, if we consider the case M = 4, R = 3 we obtain 

= ^+Pi) + ai(2ga) + a2(4 g | + g ?) + 4a 3 (2 9 | + 9 ?g2) 

, m« 4 , 10 2 2, 4x, ^ , 2g 2 , (4gi + g?) , (8g| + 4g 2 g 2 ) 
+a 4 (16g 2 + I2 qi q 2 + g ) + _ + 7l _ + 72 _ + 73 _ 

qf q{ ql q\ 

(2.22) 

where the notation in [23j corresponds to A = 2£o,7j = 2 2j+1 £j (J = 1,2,.. .). As it can be 
easily appreciated, M is just the degree of the polynomial perturbation and (2R + 2) is the 
maximum degree in the denominator of the rational perturbations. 
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3 Hamiltonian systems on the Poisson algebra sZ(2,R) hz 



The aim of this Section is to show that all the abovementioned integrable systems can be 
written within a common algebraic framework as Hamiltonian systems on the six-dimensional 
Poisson algebra s/(2,R) © /t 3 , with generators {J+, J-, J 3 , A + , A-, M} and non- vanishing 
Poisson brackets given by 

{J 3 ,J+} = 2J+ {J 3 ,J_} = -2J_ {J_,J+} = 4J 3 {A^,A + } = M (3.1) 

(all the brackets between generators of si (2, R) and /i 3 vanish). Obviously, this Poisson algebra 
has two Casimir functions: the central generator M and the sl(2, R) Casimir function 

C s i(2) = J+J- - Jl- (3.2) 

As we shall demonstrate in the next Section, this underlying symmetry of the integrable 2D 
HH systems will provide ND generalizations for all of them in a systematic way. 

In particular, let us consider the two-particle symplectic realization of sZ(2,R) © /i 3 given 

by 

J+ = Pi J- = Qi J3 = 9iPi (3.3) 
A + =p 2 A„=q 2 M = 1. (3.4) 

whose symplectic leaves will be labelled by the values of the Casimir functions, namely 
Csi(2) = an d M = 1. Now it is easy to realize that all the 2D integrable HH Hamiltonians 
previously described are of the form 

^ = l -(pi+pl)+V{qlq 2 ). (3.5) 

Now it becomes clear that is sl(2, R) © /i 3 symmetric, since it would be just the N = 2 
realization (|3.3p and (|3,4p of the generic Hamiltonian 

U = \{J + + A 2 + )+V{J-,A^). (3.6) 

Moreover, for each specific choice of T~L that corresponds to a given integrable HH system, it 
turns out that we can find an "abstract" function X = I(J+, J-, J 3 , A + , A_,M) in involution 
with ~H with respect to the Poisson bracket (|3.ip and whose symplectic realization (|3.3p ~ (|3.4p 
is just the associated integral of the motion I^ 2 \ In fact, the only thing that we have to do 
is to write the corresponding 1^ as the most generic function of the form 

and perform the substitution 

1 {2) (pI, Qi, ?i Pi, 92, P2) ^1(J+, J-, Js, A + , A_). (3.8) 

Here by 'most generic' substitution we mean that, for instance, terms of the type q 2 Pi have 
to be replaced by a linear combination of the type a\ J+ J- + a 2 J$ ■ Finally one has to impose 
that the function I so obtained does Poisson commutes with (|3.6p by using the bracket (|3.ip , 
and this process will fix unambiguously all the constants aj. In this way the 'abstract' integral 
of the motion I is obtained, thus obtaining explicitly the full s/(2,R) © /i 3 invariance of the 
system. As we shall see in the sequel, this scheme can be sucessfully applied onto all the 
integrable HH systems that have been previously described. 
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3.1 The Sawada-Kotera case 



If we consider the N = 2 realization (|3.3p and (|3.4p of the algebra sl(2, R) © /13 we can write 
the perturbed Sawada-Kotera Hamiltonian (|2.ip as 



(3.9) 



and the integral of the motion ([2.20 can be writen as an 'abstract' sZ(2,R) © /13 funtion in 
the following way: 



-J+A% + 2SA-(5A-J- + A+ J 3 ) + 2aSA-J- 



' ' 1 2 1 OX/1 f X A 7 i/l T_ \ 1 0«,Jf /I T I ^ + ^4 2 



2J3 
3 



+a aAtJ_ 



A 2 _ J_ 
1T + T 



o 



A+ A. 



+ A[2A_J + +A + J 3 ] 



-4? 



4a 



(3.10) 



By making use of the Poisson algebra (j3.ip . a straightforward computation shows that ([3.90 
and ([3.100 are in involution, and this will be true for any further symplectic realization of 
the sZ(2,R) © /13 Poisson algebra that we could consider. 



3.2 The Kaup-Kupershmidt case 

In this case the very same procedure for ([2.30 gives us the sZ(2,R) © /13 invariant object: 



H 



X - ( J+ + A\) + «5(J_ + 16Ai) + a + y A 3 ) + A + 



(3.11) 



and the integral of the motion coming from ([2.40 is proven to be 
X 



\j\ + 5(J- [35 J- + J+] + 2J 3 2 ) 



+q J_(A_ J + - A + J 3 ) - a J 2 



h A 2 

6 



+ 2A„(J 3 2 -8J 2 _ 



+A ( j- (j+ + 2aA -) + J> ( 2aA - + 25 + T 

Once again, the involutivity between ([3.120 and ([3.110 can be easily checked 



2A v 



(3.12) 



3.3 The KdV case 

In this case the HH Hamiltonian (|2.5[) is given by 

-H = \{J+ + A\) + 5(J- + A 2 _) + QA 2 _ + a(J_A_ + 2A 3 _) + A (3.13) 
and ([2.60 is just the symplectic realization of the 'abstract' integral 

X = 5(^J + + {35-n)J^-^J + + a(-A^J + +aJ-(^+A 2 _\+A + J 3 



+2a5A_J- + -^(35-0,- 2a A. 



(3.14) 



In the same way, the full perturbed KdV system (|2.20p is written in terms of the Poisson 
algebra generators as 

1 A M R 

H M ,R = 2( J + + ^ + )+l- + E«^ + Ew^ 

i=l i=l - 

M III /, , A 

tA; At— 2k 



i=l fc=0 V J 

+ EE^- 2fc ( 2 "/)^ (M>i?) (3.15) 



i=l fc=0 

and the integral of the motion coming from (|2.21|) is found to be 

2A 

1m,r = -A.J+ + J 3 A+- —A- 



M 



jk *i—l—2k 



\ i=l fc=0 V ' 

4 A^-dimensional symplectic realizations 

The previous sZ(2, R) ffi/13 Hamiltonian structure of the HH systems will allow us to construct 
ND integrable generalizations of them in a systematic way. In particular, we can firstly 
introduce the following ND symplectic realization of s/(2,R) © /13, in which the first set of 
(N — 1) degrees of freedom is associated to sl(2, R), while the last one is just the usual 
one-particle I13 realization with M = 1: 

N-l N-l N-l 

J + = E? 2 J - = E? 2 ^ = 52^ 

1=1 1=1 1=1 

A + =p N A_=q N M = 1. (4.2) 

It is straightforward to prove that any ND Hamiltonian constructed as the previous symplectic 
realization of any function Ti on sZ(2,R) © /13, namely 

%W =n(J + ,J-,J 3 ,A + ,A_,M) (4.3) 

is quasi-integrable, since the following (N — 2) functions 

m 

C M = J2 (q iPj - q jPl f m = 2, . . . , N — 1 (4.4) 

l<i<j 

are constants of the motion for Ti, and they are mutually in involution. In fact (|4.4p are just 
the so-called sl(2, R)-coalgebra integrals for the sl(2, R)-sector of the Hamiltonian (see [251 
26, 27, 28j for details). Note that all these constants of the motion depend on the first (N — 1) 
canonical variables, and any contribution in the Hamiltonian coming from the last degree of 
freedom (qN,PN) will obviously Poisson-commute with any of them. 
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Therefore, any system given by (|4.3p is only 'one integral away' from being completely in- 
tegrable. Moreover, if for a certain choice of T~L we are able to find one additional independent 
integral of the 'abstract' form 

I = I(J+,J-,J 3 ,A+,A-,M) (4.5) 

and we realize it under the ND symplectic realization (|4.ip and (|4.2p . the complete integra- 
bility of T-L^ will be guaranteed in the ND case. This will be the case for all the HH systems, 
whose 2D integrability will provide such function I. 

Thus, the ND analogue of a given 2D HH system will be defined as the iV-degrees of 
freedom realization (I4.1|) - (l4.2p of the corresponding s/(2,M) © /13 Hamiltonian (13.61) . and the 
TV-particle realization of the additional integral X will provide the complete integrability of 
the system. 



4.1 A KdV example 

In order to illustrate this general result, let us consider (I3.13P as the algebraic definition of 
a ND HH-KdV Hamiltonian. Namely, by taking N = 3 in (|4,ip - (l4.2p and, afterwards, by 
substituting these expressions into (|3.13p we get: 



J =1 



Q3 + 2gf + 



(4-6) 

Note that in this way we are obtaining a 'radial' generalization of the KdV potential in the 
{0.1^0.2) subspace, which is a direct consequence of the spherical symmetry of the symplectic 
realization (|4.ip for the sZ(2,R) sector. 

In general, the specialization of (|3. 13|) and (|3. 14[) through the iV-particle symplectic re- 
alization (|4.1|) ~ (|4.2p will provide the explicit form of an integrable ND generalization of the 
KdV Hamiltonian, that reads 



'N-1 



N-1 



o 



N-1 

E« 

i=i 



q N + 2q N + ■ 



A 



(4-7) 

and the explicit form of the integral of the motion coming from (|3.14p is immediately obtained: 



N-i 



N-1 

3=1 



-a 



N-1 



N-1 

*v Z) p2 3 + a [ Z i 2 j 



3=1 

'n-i 



3=1 





n 


)- 


2 






£ 






4 



N-1 



3=1 



+PN 



'N-1 
.5=1 



+2a5q N Qj + 



A 



'7V-1 x 



-(35 - n - 2a<?Ar) 



(4.8) 
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It is straightforward to check that (|4.8|) Poisson-commutes with the full set of si (2, R) integrals 
QUID and, obviously, with ([47f]) . 



4.2 Adding centrifugal terms 

But a further generalization of this construction is still possible by considering a more general 
symplectic realization of the Poisson algebra s/(2,R) © /13. Such generalization breaks the 
radial symmetry in the s£(2,R) sector by adding centrifugal terms, albeit by preserving the 
complete integrability of the new system (note that, in general, this is by no means guaranteed 
when centrifugal/monopole terms are added Namely, let us consider the functions 

N-i , 7 \ N-l N-l 

J + = E(^ + 4) j - = E^ 2 J s = ^2qiPi (4-9) 
i=i v 1% J i=1 i=1 

yl + =PAr ^4- = g7V M = l. (4.10) 

which define the (&i, £>2, • • • > fr/v-i) © (1) symplectic realization of s/(2, R) 0/13. Note that its 
only difference with respect to (|4.ip and (|4.2p is the addition of (iV — 1) centrifugal terms in 
the realization of the J + generator (despite the fact that the constants bi can be negative, we 
shall continue calling those terms as 'centrifugal' ones). 

When these terms are included, the constants of the motion (|4.4p coming from the si (2, R)- 
sector of the Hamiltonian are converted into (26^ [28j [30] 

m m / 2 2\ 

C M= (m-QjPif+ E 4 + l ''7 m = 2,...,iV-l. (4.11) 

l<t<j l<i<j V ^ ^ / 

But it is essential to stress that any sZ(2,R)©/i3 Hamiltonian ()4.3p . when written in terms of 
the symplectic realization (|4.9p - (|4.10p will Poisson-commute with all the integrals (|4.11|) . And 
the same will happen to the corresponding symplectic realization of the additional integral 
I, that will be modified only through the centrifugal terms coming from the J+ generator. 

Namely, if we consider the N = 3 generalized symplectic realization (|4.9|) — ()4. 10p of the 
KdV Hamiltonian (I3.13p . we obtain 




q3 + 2ql I +- 



A 



=1 

(4.12) 

This system is, by construction, completely integrable: its integrals of the motion in involu- 
tion are the integral given by (I4.1ip and the N = 3 generalized symplectic realization 
(|4.9p - (|4.10p of the 'abstract' integral (|3.14p . Note that this Hamiltonian has two types of 
'centrifugal' terms, one of them linked to the parameter A and the other one induced by 
the generalized symplectic realization and controlled by the parameters 6j. We recall that 
when A = the ND generalization of (|4.12p has been shown to be integrable even when the 
harmonic oscillator term is fully anisotropic (see [19] and references therein). However, in 
our case the sl(2, R) /13 symmetry imposes the isotropy in the first (N — 1) coordinates, 
thus allowing the existence of the 'global' centrifugal term linked to A. 



11 



5 ND integrable Henon-Heiles systems 



Finally, we present the complete expressions for the most general HH systems that we can 
obtain through the algebraic procedure described in the previous Sections. We stress that 
the following three ND Hamiltonians do commute with the same 'universal' set of integrals 
(|4.1ip . Also, it becomes evident that finding in each case the remaining integral l( N > without 
the aid of the underlying sl(2, R) © /13 symmetry would be a quite complicated task. On the 
other hand, this symmetry guarantees by construction that the integral X^ N ' is in involution 
with the set of functions (14.11i . 



5.1 The generalized ND system of Sawada-Kotera type 



When the symplectic realization (|4.9p ~ (|4.10p is applied onto (|3.9p . the ND HH-Sawada-Kotera 
Hamiltonian reads 



N-l 



« <w) = 5 E 



i=l 



Z 1 1 

Pi + -3. 



'N-l 



a qN 



i=l 



'N-l 
. i=l 



+ Iqn I + 



A 



and the integral of the motion coming from (|3.10p now reads: 



N-l 



i=l 



N-l 



N-l 



Pi 



i=l 



+2aSq 



N 



N-l ' 

E* 2 

. 1=1 

/N-l 



N-l 



QN 



+a \aq 2 N IJ2 q\ 



1=1 

a 2 1 ^ 



i=l 



AT-1 



/JV-1 



+f IE« 



iV-l N 
2 



i=l 
AT-1 

E 

8=1 



i=l 



'AT-1 



+ ^(E 



8=1 



1=1 



2^AT ^2 \Pi + ~\ ) +PN ( E 

i=l V 9i / \ - =1 



+A 



4a 



P% , 
. i=l 



/ 



/AT-1 N 

V t=i / 



(5.1) 



(5.2) 



The fact that both quantities are in involution can be straightforwardly checked by direct 
computation. 
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5.2 The generalized ND system of Kaup-Kupershmidt type 

In the same way we can obtain the generalized version of this system by considering the 
symplectic realization (|4.9p - (|4.10p of the 'abstract' Kaup-Kupershmidt Hamiltonian (|3.1ip 



'N-l 

E 

V 1=1 



+a q N 



Pi+ — 
if 

N-l ' 



'N-l 



i=i 



E 

i=i 



16 



+ Y qN r 



tz 1 \ /N-l \ ' 



(5.3) 



In this case the integral of the motion (I3,12|) reads 



N-l 



!(E 



. i=i 
'n-i 



1=1 



Z i * 



'AT-1 



N-l 



N-l N-l 



+« [ E q i E (y p2 i + ^ J ~ PN E qip 

N-l \ 2 /AT-1 \ 



3*E^ 2 + E(^ 



i=l 
AT-1 



i=l 



i=l 



'N-l 



+ 2 E 



1=1 



a 



E- 

\ i=i 



JV-l 



1 \ ^ 2 2 

r 2^ ?i + ?jv 



i=l 



+2ag/v 
/ 



+A 



/N-l 

\i=i 



JV-l 



i=l 



E(P? + 3) + 



A 



+- 



3u 



/N-l 
\i=l 



2agAr + 25 + 



1 



JV-l 

i=i 



N-l 



\ 



/N-l 

( E 

V i=i 



i=l 



+ 



2A 



AT-1 

E 9. 
1=1 



+ 



/AT-1 

E^ 2 

V i=i 



5.3 The generalized ND KdV system 



Finally, the generalized version of the T~i^ R Hamiltonian is obtained by substituting the 

lization : 



generalized symplectic realization into f)3. 15|) . The final expression is 



n 



(TV) 
M,R 



P 3 ^ „2 



+ Pn + 



A 



N-l 



M 



i=l k=0 
R If] 

+EE^~ 



-2k — h 
k 



'N-l 

E^l 

.3=1 



i-2k 



2k 



i=l k=0 



i — k 



In 



2k 



'N-l " 



i+l-k 



(M > R) 



(5.5) 
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and the same procedure for (|3.16|) gives rise to the additional integral 



r (N) 

-M,R 



N-l 



N-l 



-in E ( P 2 j + ^ J + pn [ E Vi<h 



2A 



JV-l 



rt=ll 
2 J 



/ M [ 

+ (E^ 2 I EE«» 2 



i-\-2k 



JV-l 

E 



-<?/v 



Ov-1 



2=1 fc=Q 



r m 



E E 



2A- 



i=l fc=0 



i — l — k 
k 



i + l-k 



'N-l 



E^ 2 I ti 1 

3=1 



-2fc 



■i+l-2fc 

In 



V 



'AT-1 



i+2-fc 



E*l 



(M > 



(5.6) 



6 Further generalizations 

Summarizing, the main result contained in the previous pages is the fact that the multipara- 
metric ND Hamiltonian 

nW = \ E (*? + |) + \pn + V (E fl?, ^ (6.1) 

is completely integrable provided that the function V is any of the 2D integrable HH poten- 
tials. Moreover, the full set of integrals of the motion for (|6.ip can be explicitly obtained as 
a consequence of its underlying sl(2, R) © /13 symmetry. 

However, none of the arguments underlying the previous construction prevents us to use it 
for any 2D integrable potential of the form V (x 2 , y) (see the classifications in [T5 | l20 j l33" l [34] ) . 
For instance, the 2D Holt potential (see [201 [34]) 

n^ = \{pl+pl) + q ^ + §d) (6.2) 

is also of the form V (x 2 , y) and its (quartic) integral of the motion is given by 

X = pf + 2p 2 1 p 2 2 + 24ql P2 q lPl +4qPq 2 1 pl + 72qlql (6.3) 
Therefore, this system admits a sZ(2,K.) © formulation, that leads to the Hamiltonian 

U = \{J + + A%) + A~J (j_ + 9 -A 2 _}j (6.4) 

and to the abstract integral 

X = Jl + 2J + A 2 _ + 2AAlA + Jz + 4:A~JJ + J- + nA^J-. (6.5) 

As a consequence, by substituting (|4.9|) - (|4.10p in (|6.4|) and (|6.5p we get an ND integrable 
generalization of the Holt potential, for which the remaining integrals are again given by 
the 'universal' set (|4.11|) . As a side remark, in case that the additional integral X would be 
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unknown, note that this algebraic approach allows the efficient use of symbolic computation 
packages in order to look for such a constant of the motion, since the ND Hamiltonian can 
always be written as an object defined on the set of six generators that generate the Poisson 
algebra s£(2,R) © h 3 . 

Besides its immediate applicability to many other Hamiltonians, the Poisson symmetry 
framework here presented opens the path to several further developments on HH systems 
based on such algebraic perspective, that we plan to work out in the near future. Firstly, 
the quantization of all the integrable HH systems here presented can be faced by considering 
the Lie algebra s/(2,R) ffi/13 as the corresponding underlying symmetry (regarding the many 
facets of the quantization problem for the original HH system (II. ip see [3]). Therefore, the 
explicit solution of these quantum HH systems should be amenable by making use of the 
representation theory of this algebra and by diagonalizing the Hamiltonian % in a common 
eigenbasis with respect to the (suitable defined) quantum 'abstract' invariant X. 

Secondly, the type of Poisson approach here used should be the suitable one in order to 
define appropriate (integrable) analogues of the HH systems on curved spaces, that should 
lead to the flat HH Hamiltonians when the curvature parameter is turned off. In this respect, 
we have recently presented in [30J a constructive si (2, R)-Poisson algebra approach to many 
integrable systems defined on ND conformally flat spaces. In particular, the problem of 
getting a consistent definition of integrable HH systems on the spherical and hyperbolic 
spaces should be guided by the preservation of some underlying symmetries that should be 
related to the flat sl(2, R) © /13 ones, and constitutes a challenging problem. 

Also, it is well known that one of the main features of the Hamiltonian systems endowed 
with a given Poisson coalgebra symmetry is the fact that integrable deformations of them 
can be systematically constructed by introducing a suitable (/-deformation of the underlying 
coalgebra (see |25| . I26| . l28j ^ . In our case, one could assume the same 'abstract' form of the HH 
Hamiltonians here considered but -for instance- with the functions (J+, J_, J3) being now 
generators of a given (/-deformation of s/(2,R). In that case, the quantum group theory tells 
us how to deform the full set of integrals (|4.1ip . and the only remaining task would be to find 
the appropriate deformation of the integral I. Previous experience on the subject shows that 
integrable systems on curved spaces can arise in this way, and in such cases the deformation 
parameter is directly related with the curvature of the underlying space |38| . 

Finally, we recall that (1+1) Minkowskian HH systems with kinetic energy of the type 
ip\ ~ P2) have been considered in the literature in connection with geodesies in gravitational 
waves [39]. Indeed, the integrability properties of this class of HH systems could be anal- 
ysed in terms of the sl(2, R) © /13 symmetry, thus allowing their generalization to an (N+l) 
Minkowskian space by following the approach here presented. Another interesting problem 
in a General Relativity context seems to be the study of the algebraic integrability properties 
of HH systems perturbed by a Kepler potential [3D], a term that in our language would be 
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written as a perturbation of the type + A 2 _. 




15 



References 

[1] M. Henon, C. Heiles, Astron. J. 69, 73 (1964). 

[2] M. Tabor, Chaos and integrability in nonlinear dynamics, (Wiley: new York) (1989). 

[3] M.C. Gutzwiller, Chaos in Classical and Quantum Mechanics, Interdisciplinary Applied 
Mathematics, Vol. 1 (Springer: New York) (1990). 

[4] D. Boccaletti, G. Pucacco, Theory of Orbits, (Berlin: Springer) (2004). 

[5] T. Bountis, H. Segur, F. Vivaldi, Phys. Rev. A 25, 1257 (1982). 

[6] Y.F. Chang, M. Tabor, J. Weiss, J. Math. Phys. 23, 531 (1982). 

[7] B. Grammaticos, B. Dorizzi, R. Padjen, Phys. Lett. A 89, 111 (1982). 

[8] J. Hietarinta, Phys. Rev. A 28, 3670 (1983). 

[9] A.P Fordy, Phys. Lett. A 97, 21 (1983). 

[10] S. Wojciechowski, Phys. Lett. A 100, 277 (1984). 

[11] R. Sahadevan, M. Lakshmanan, J. Phys. A: Math. Gen. 19, L949 (1986). 

[12] A.P. Fordy, Physica D 52, 204 (1990). 

[13] W. Sarlet, J. Phys. A: Math. Gen. 24, 5245 (1991). 

[14] V. Ravoson, L. Gavrilov, R. Gabor, J. Math. Phys. 34, 2385 (1993). 

[15] M. Lakshmanan, R. Sahadevan, Phys. Rep. 224, 1 & 2 (1993). 

[16] R. Conte, A.P. Fordy, A. Pickering, Physica D 69, 33 (1993). 

[17] M. Salerno, V.Z. Enol'skii, D.V. Leykin, Phys. Rev. E 49, 5897 (1994). 

[18] C. Verhoeven, M. Musette, R. Conte, J. Math. Phys. 43, 1906 (2002). 

[19] N. A. Rostov, V.S. Gerdjikov, V. Mioc, J. Math. Phys. 51, 022702 (2010) 

[20] B. Grammaticos, B. Dorizzi, A. Ramani, J. Hietarinta, Phys. Lett. A 109, 81 (1985). 

[21] M. Antonowicz, S. Rauch-Wojciechowski, Phys. Lett. A 163, 167 (1992). 

[22] A.N.W. Hone, V. Novikov, C. Verhoeven, Inv. Probl. 22, 2001 (2006). 

[23] A.N.W. Hone, V. Novikov, C. Verhoeven, Phys. Lett. A 372, 1440 (2008). 

[24] B. Grammaticos, A. Ramani, Phys. Lett. A 139, 299 (1989). 

[25] A. Ballesteros, O. Ragnisco, J. Phys. A: Math. Gen. 31, 3791 (1998). 

[26] A. Ballesteros, F.J. Herranz, F. Musso, O. Ragnisco. Superintegrability in Classical and 
Quantum Systems CRM Proceedings and Lecture Notes 37 ed P Tempesta et al (Prov- 
idence, RI: AMS), 1 (2004). 

[27] A. Ballesteros, A. Blasco, J. Phys. A: Math. Theor. 41, 304028 (2008). 



16 



[28] A. Ballesteros, A. Blasco, F.J. Herranz, F. Musso, O. Ragnisco, J. Phys: Conf. Ser. 175, 
012004 (2009). 

[29] A. Blasco, Integrability of non-linear Hamiltonian systems with N degrees of freedom, 
PhD. Thesis, Burgos University (2009) 

[30] A. Ballesteros, A. Enciso, F.J. Herranz, O. Ragnisco, Ann. Phys. 324, 1219 (2009). 

[31] R. Conte, M. Musette, C. Verhoeven, Theor. Math. Phys., 144, 888 (2005). 

[32] B. Grammaticos, B. Dorizzi, A. Ramani, J. Math. Phys. 25, 12 (1984). 

[33] J. Hietarinta, Phys. Rep. 147, 87 (1987). 

[34] A.M. Perelomov. Integrable systems of classical mechanics and Lie algebras, (Basel: 
BirhauserVerlag) (1990). 

[35] G. Tondo, J. Phys. A: Math. Gen. 28, 5097 (1995). 

[36] A. Ramani, B. Dorizzi, B. Grammaticos, Phys. Rev. Lett. 49, 1539 (1982). 

[37] A. Ballesteros, F.J. Herranz, J. Phys. A: Math. Theor. 42, 245203 (2009). 

[38] A. Ballesteros, A. Enciso, F.J. Herranz, O. Ragnisco, Phys. Lett. B 652, 376 (2007). 

[39] K. Vesely, J. Podolsky, Phys. Lett. A 271, 368 (2000). 

[40] F. Kokubun, Phys. Lett. A 245, 358 (1998). 



17 



